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TECHNICAL NOTE 3 ^ 3 ^ 

ON THE KERNEL FUNCTION OF THE INTEGRAL EQUATION 
REIATING LIFT AND DOWNWASH DISTRIBUTIONS OF 
OSCILLATING WINGS IN SUPERSONIC PLOW 
By Charles E. Watkins and Julian H. Bennan 


SUMMARY 


This paper treats the kernel function of the integral equation that 
relates a known or prescribed downwash distribution to an unknown lift 
distribution for harmonically oscillating wings in supersonic flow. The 
treatment is essentially an extension to supersonic flow of the treatment 
given in NACA TN 3131 for subsonic flow. For the supersonic case the 
kernel function is derived by xise of a suitable form of acoustic doublet 
potential which employs a cutoff or Heaviside unit function. The kernel 
functions are reduced to forms that can be accurately evaluated by con- 
sidering the functions in two parts: a part in which the singularities 

are Isolated and analytically expressed, and a nonsingtOar part which 
can be tabulated. 

The kernel is treated for the two-dimensional case, and it is shown 
that the two-dimensional kernel leads to known lift distributions for 
both steady and oscillating two-dimensional wings. The kernel function 
for three-dimensional supersonic flow is reduced to the sonic case and 
is shown to agree with results obtained for the sonic case in NACA 
TN 3 I 5 I-# aJ3d the downwash functions associated with "horseshoe" vortices 
in supersonic flow are discussed and ejqjresslons are derived. 


INTRODUCTION 


In reference 1 the kernel function of an Integral equation relating 
a known or prescribed downwash distribution to an imknown lift distri- 
bution for a harmonically oscillating finite wing of arbitrary plan form 
was treated for compressible subsonic flow. The purpose of the present 
paper is to extend this treatment of the kernel function to supersonic 
flow. 


Qlie kernel functions under consideration arise when llnearlzed- 
boundary-value problems for obtalnlxig aerodynamic forces on oscillating 
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vlngs are reduced to Integral equations involving the distribution of 
pressure or wing loading as the iinknown. In such integral equations 
the kernel functions play the important role of aerodynamic Influence 
functions in that they give the normal induced velocity or downwash at 
any one point in the plane of the wing due to a unit pressure loading 
at any other point in the plane of the wing. 

As the kernel fxmctlons arise in the analysis, they eire mathemati- 
cally defined by rather intricate improper integrals and possess singu- 
larities as high as second order. It is therefore desirable to Isolate 
the singularities and deteimlne their explicit nature in order to make 
the integral equation more amenable to solution, in particular amenable 
to solution by approximate or numerical procedures. 

Approximate lifting-surface theories for finite wings, such as the 
methods developed by Falkner and Multhopp (refs. 2 and 5) and others, 
have afforded considerable success in the calculation of aerodynamic 
coefficients for steady subsonic aerodynamics. Similar approximate 
methods have been successfxilly employed to obtain coefficients for 
two-dimensional oscillating wings in subsonic (compressible) flow (for 
example, refs. 4 and 5) and are now being extended to the finite oscil- 
lating wing in subsonic flow by Harry L. Runyan and Donald S. Woolston 
of the langley Aeronautical laboratory and by ¥. P. Jones (ref. 6) . It 
is reasonable to expect that these methods can be further extended to 
apply to finite wings in supersonic flow. 

In supersonic flow, solutions of the boundary- value problem for 
some particular plan forms and downwash conditions can be obtained in 
the form of Infinite series in terms of a parameter involving the fre- 
quency of oscillation (see, for exan^jle, refs. 7 "to lO) or in the form 
of rather complicated definite Integrals (refs. 11 and 12). The infinite- 
series method furnishes a relatively simple means of obtaining the loading 
on oscillating wings for low values of the frequency parameter, but for 
large values of this parameter the series expansions converge so slowly 
that recourse must be had to other proced^n:eB for obtaining the wing 
loading. One feasible method is to study and develop approximate pro- 
cedures for solving the Integral equations that Involve the unknown 
loading and its associated kernel function. The first step toward such 
a development is to isolate and determine the explicit natinre of the 
singularities of the kernel fimctlon; this step is accomplished in the 
present report. 

The report contains the derivation of the kernel function in the 
form of an improper integral and a reduction of this Integral to proper 
form. The singularities of the kernel function are Isolated and expressed 
analytically, and the nonsingular parts are reduced to a form readily 
amenable to nunterlcal evaluation, as was done in reference 1 for subsonic 
flow. Seme expanded forms of the kernel function are derived, and one of 
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these is izsed to ohtaln a reduction to two-dimensional flow. In 
appendix A, the limiting case for sonic flow is derived and shown to 
agree with the res^llts in reference 1. Appendix B is devoted to certain 
integrals of the kernel function. Oliese integrals relate to "horseshoe" 
vortices in supersonic flow, as treated, for example, in the steady case 
hy Schllchtlng in reference I 3 , and may he of Interest in certain modes 
of application. 


SYMBOIS 


Il(x) 

Jn(x) 

Kl 

K(Xo,yo) 

k 

K (xo) 

L(|) 

Ll 

I 

M 

P 


velocity of sound 

Bessel function of first kind, of first order with 
Imagineiry argument 

Bessel function of first kind, of order n 

modified Bessel function of second kind, of first order 

kernel function for three-dimensional flow 

reduced-frequency parameter, iq/V 

kernel function for two-dimensional flow 

lift dlstrlhutlons 

modified Struve function of first order 
unit length 
Mach number , V/c 
perturbation pressure 


r 


]Jyo^ + 


s 

t 

U(x) 


region of xy-plane occiq)ied by wing 
time 

unit function 


V 


forward velocity of wing 
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w(x,y,t) 

w(x,y) 


x,y,z 


Xq = X - I 


downwash velocity, e^^(x,y) 

coc^lex araplitiode function of prescribed vertical 
velocity 

Cartesian coordinates attached to wing moving In 
negative x-dlrectlon 


yo = y “ 
p = - 1 

S(x) 

P 

JZ5(x,y,z,t) 

0(x,y,z) 

t(x,y,z,t) 

if(x,y^z) 

cu 

CO = co/V^2 


Dirac delta fimctlon 

Cartesian coordinates used to represent space location 
of doublets 

fl\ild density 

velocity potential, e^‘^^(x,y,z) 

congjlex angjlltude function of velocity potential 

acceleration potential, eiti>ti^(x,y, z) 

complex amplitude fimctlon of acceleration potential 

clrcxilar frequency of oscillation 


ANALYSIS 

Integral Equation Eelatlng Downwash and Lift Distribution 


Ihe Unearlzed-boundary-value problem for the determination of the 
aerodynamic forces on a wing can be Immediately reduced to a problem of 
solving an Integral equation that relates downwash and lift distribution. 
Hhe purpose of this section Is to Introduce and briefly discuss this 
equation. 

Since the Integral equation has the same formal appearance for sub- 
sonic and supersonic flow and Is derived In various publications (for 
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example, refs. 1 btk^ llj.), the equation will not he rederived here hut 
will he formally stated so as to serve as a starting point in the 
analysis. In keeping with linear theory, the wing is considered as a 
plane^ in 5 >enetrahle surface S which lies nearly in the xy-plane as 
indicated in the. following sketch: 



TTip x,y, z coordinate system and the siarface S are assumed to move in 
the negative x-direction at a -uniform veloci-ty V. 

In tenns of these coordinates, the integral eq.ua tion may he formally 
written as 

w(x,y) = JJ L(|,Ti)K:(x-|,y-Ti)d| dn (l) 

where w(x,y) is the complex amplitude fimction of prescribed vertical 
velocity or downwash at points in S and is defined as follows; 

w(x,y,t) = e^’^(x,y) 

where cu is the frequency of pulsation or oscillation. !Ehe kernel 
fimction K(xQ,yo) physically represents the contrihutlon to the down- 
wash at a point (x,y) in S dxie to the presence of a pulsating pressure 
doublet of unit strength located at some other point (|,ti) in S. It is 
a fxmction not only of x, y, and q, hut also of Mach number and 
frequency. The function L(|,ti) in eqiiatlon (l) is the vmknown lift 
distribution or local doublet strength. (Althovi^ it is usually conven- 
ient to factor out the density term l/p as indicated in eq. (1), 
this was not done in ref. 1.) 

Equation (l) pertains formally to either subsonic or supersonic 
flowj however, separate treatments of the two cases are reqiilred because 
of -wide differences associated with flow characteristics. So far as the 
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Integral, equation is concerned, the differences in the two cases lie 
mainly in the ]sernel functions; These differences are associated with 
the differences in character of doublets for the two cases. Although 
the main purpose of this analysis is to derive and treat the kernel 
function for the supersonic case, a necessary first step is to formulate 
a doublet s\ii table for such a treatment. In the following section, a 
desired form, which was arrived at by a convenient use of a cutoff, or 
Heaviside unit function, is presented in equation (5) • 


Pulsating Doublet Moving at Svipersonic Speed 

The governing differential equation for linearized unsteady flow 
at either subsonic or supersonic speeds, which the doublet potentials 
must satisfy, is the well-known wave equation referred to a moving 
coordinate system: 

^ ^ ^ _ -L/V = 0 

8y^ 8z^ c^V htJ 

Under the assumption that dlstiurbances vary harmonically with respect 
to time, this equation becomes 



where is a coorplex amplitude function defined by 

i|f(x,y,z,t) = e^t(x,y,z) (5) 

It may appear that, since the same differential equation (eq. ( 2 )) is 
Involved, a logical way of obtaining the potential for a pulsating 
doublet moving at supersonic speed is by simple analogy or continuation 
from the potential for the doublet moving at subsonic speed. This pro- 
cedvire is applicable only in a broad sense because, as discussed in 
reference I5 with regard to sources in supersonic flow, the potential 
of a doublet moving at supersonic speed consists of the sum of two 
effects corresponding to a retarded-type potential and an advanced-type 
potential which relate to the two wave fronts encountered by a point at 
any time; whereas for subsonic speed only the retarded type of potential 
is admissible. (The advanced- type potential for subsonic or sonic speed 
does not satisfy the Somraerfeld radiation condition, which requires that 
disturbances be propagated away from their point of origin.) In the 
second place, the potential that may be obtained by analogy with the 
potential for subsonic speed must, as subsequently discussed, be rather 
severely restricted before it mathematically describes the physical 
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realities of a disturbance moving at supersonic speed.. In the following 
development, a desired form of the doublet potential is arrived at by- 
consideration of these restrlcticsns applied to bo-th a retarded and an 
ad-vanced "type of po-fcential that may be obtained by analogy wl-th resul-bs 
for subsonic flow. 

By analogy -wi-th results for subsonic speed (for example, eq,. (A 9 ) 
of ref. 1) or, more directly, from "the discussion of source potentials 
in supersonic flow (ref. 15), the s-um \lfj3 of the re-barded and ad-vanced 
tsrpes of po-bentials required to form -the doublet potential for supersonic 
speeds may be -written wi-fch -fche doublet situa-ted at -the origin as 

V cp2 cp2 / _ \ cp2 cp2 

■ + 7============: — 

/x^ - 3^y^ - /x2 _ p2y2 _ p2^2 




gi(art;-I^fflx)QQg^^l^2 _ p 2 y 2 _ p2^2 ^ 

sj-x^ ~ 


( 4 ) 


where M = V/c, p = - i, and S = a>/vp^. The restrictions that must 

be placed on -this expression are: (a) only real -values of the radical 

-term \Jx^ - 3%^ - p^z^ sire to be considered and (b) -the -values of the 
expression and its derivatives are to be considered zero when x is 
negative. These restrictions follow from -the physical consideration 
that small disturbances propaga-te at sonic speed, and in a sxqpersonic 
stream do not progress forward of -their point of origin. 


A convenient way of writing the expression for -<(^33 -with these 
restrictions accounted for, as pre-viously mentioned, is to employ a 
cutoff or xmit function as a factor. Th-us, if represents -the 

restrlc-ted -value of ilri), -the an^ilitude f\mction of ijfQ may be -written 
as 


= 2 cos( ^fe -_p 3£ -_p2z 2) 

|/x2 - p2y2 - p2z2 


( 5 ) 


where 


U^x-p\/y^+z2 j 
U^x-pjy^+z^) = 0 


(x > pjy2 + z^) 

( X •< p^y^ + z 2 j 


( 6 ) 
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and only positive values of the radical + z2 aj-g considered. A 

method whereby this form of potential can he determined in a more direct 
manner is discussed in detail in reference l6. The discussion in this 
reference is in connection with the Green functions associated with the 
dispersion of sound waves in an n-dlmensional mediian in which a pulsating 
source exists. When appropriate changes are made in notation, the 
expression for given in equation (5) agrees essentially with results 

for the dispersion of waves in a three-dimensional space given in equa- 
tion ( 55 ) > chapter XVT of reference I 6 . 

With regard to the unit function U(x), in many applications where 
this function is enqjloyed it need not be defined as having any particular 
value when its argioment is zero. 3 ji other applications, especially where 
the \mit function is involved in a Fourier analysis, it must be defined 
as having a valiae of 1/2 when its argument is zero. In the present case, 
it is conveniently defined, as may be noted in equatlrai (6), as having 
zero value when its argument is zero. 

Derivatives of the unit function give rise to an Impulse function 
called the Dirac delta function. For example, 

^ U(x) = 5(x) =0 (x / 0) 

ox 

U(x) = &(x) =00 (x = 0 ) 
ox 

A useful integral property of this delta function is 

b 

f(x) 5(x) dx = f(0) (a > 0 ) 

.a 

The next step in the analysis is to make use of the doublet 
potential (eq. ( 5 )) to derive the kernel function for supersonic speed. 


Derivation and Reduction of Kernel Function 

In this section the kernel function is derived and presented. The 
function is given in terms of an improper integral by eqmtlon (I 5 ) and 
in a reduced form with no improper integrals by equation (15). As it 
is frequently desirable to present resTilts in terms of nondlmenslonal 
length variables, the restilts given in equation ( 15 ) are presented in 
this manner in eqiiation ( 16 ) . 
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In order to derive the kemel function, the f\anction if of 
equation ( 5 ) is considered as the con^ilex amplitude of the acceleration 
potential. As such, ijr is directly proportional to a perturbation 

pressirre field p = e^^p, throu^ the simple relation 

P = -pf (?) 

and to a velocity potential 

0 = 

throxi^ the equation 

V ^ + iffi^ = ^ (8) 

ox 

differentiation of equation (8) wiidi respect to z and integration 
of the result with respect to x, the vertical velocity associated with 
the acceleration potential ijr is obtained. Thus, when ijr is considered 
as the potential of a presstife doublet, eq'uatlon (8) affords a straight- 
forward means for obtaining an eqmtion for K(xo,yo) t namely; 

K(Xo,yo) “ ^ ^(Xo#yo>2) z=0 (9) 

Details of the procedure are as follows: 

The result of the differentiation of eqmtlon (8) with respect to 
z may be written as 


8x 8z 8z 8z 


( 10 ) 


When this equation is considered as 6in ordinary differential equation 

with dependent variable ^ and Independent variable x, a con^jlete 

8z 


solution is 


^ = ie 

8z V 


icux 

V 




y^z) 


itoA 

j y dx 


(11) 


where the lower limit of integration is employed in place of a constant 
of integration and may be considered as a condition that ^ vanish far 
ahead of the origin X = 0. Bius, frcan eqiiatlon ( 9 ) there is obtained 
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pXo lojA 

K(Xo,yo) = 11m ie" V / ^ %Ck,yo>^) e ^ dA (12) 

z ^0 ^ ^-00 oz 

or, after substitution of the expression for (eq. ( 5 )) into equa- 

tion ( 12 ), the results may he written as 


K(xo,yo) = 


^ z — >0 


Icuxq 

“ V 



U(A-Pr)cos(wa^A2 - p2r2 ) 

^pr &z2 

\^A2 - p2r2 




C 15 ) 


where r = ^y©^ + and, since the Integrand is zero for X < ^r, the 
lower limit of integration has been changed from to pr. 


It is apparent upon examination of equation ( 15 ) that, if the indi- 
cated differentiation ^inder the integral sign is carried out, the 
Integrand has singular and perhaps troublesome terms. Ihe Indicated 
differentiation with respect to z, however, can be replaced by equiva- 
lent operations and followed by integrations by parts that lead to a 
reduced form of the kernel function containing no Improper integrals. 
These steps follow. 


Reduced form of kernel function .- As may be directly verified, the 
indicated differentiation with respect to z in equation (15) is, in 
the limit z — >0, identical with 


-lu5A 


lim 
z — >0 


3z2 


U( A- pr) cos ( - p^r^ ) 

\/a 2 - p2r2 


1 ^ 

,2 8a1 


yo 

a-lSA 
yo^ L 


|u(vplyo| ) 


.-ISA 


./To M 


^A2 - p2y^2 


u(A-pjyoPsin(M5/A2 - p2yo2 ) _ 1 s(A-P jyo| )Bln(McEn/A2 - P^yo^ ) - 


-■I-— . 5(A-p|yo|)cos(Ms/A2 - P^q^ ) 

|/a 2 - p2yo2 


(Uf) 
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Since the coefficients of 6 ^A-P|yo| ) in eq\aation (l^l-) vnnish at 
A = Pjyoj^ it follows from the integral properties of &(A-p|yo| ) that, 

when equation (l 4 ) is suhstltuted into equation (l 5 )j the Integrals 
involving the delta function vanish and equation (15) heccmes: 

xpe-t^o 

\Jxq^ - 

i u(xo-p|yo|)sin(j€>|/xo2 _ p2y^2 J + 

r ° u(A-p[yo|)sin(l©\fc^^~I~p^^)dA 

^ ^P|yo| 

Eq\:uatlon (15) provides an expression for the kernel function that 
Involves no inproper Integrals. Except for the Integral, the terms of 
the expression can he qiilte easily evaluated with the aid of trigono- 
metric tables, except at yg = 0 , where the function is singular, and 
at Xg = p|yo|, where the fmction is Indeterminate. Bae integral is 
well behaved and can be accurately evaluated by numerical or approxi m ate 
procedures. The singularities and indeterminate valties are isolated ^d 
discussed in a later section, but it is desirable first to express the 
function K(xg,yg) in terms of nondimensi anal length variables. As a 
check; on the correctness of equation (15)^ the expression for K(xg,yo) 
is also reduced to the Umitlng value for M = 1 and compared in 
appendix A with the corresponding limiting value for the subsonic case. 

Ihe kernel fiinction in terms of nondimensional length variables . - 
Althou^ the preceding results contain dimensional length variables, it 
is usually desirable to have such results in terms of nondimensi onal 
length variables. Ey employing the variables Xg and yg in a new 

sense to mean that they have been refeirred to some chosen length l and 
by introducing the reduced-frequency parameter k = Zcs/V, the length 
variables may be made nondlmenslonal. (in flutter theory the reference 
length normally is selected as a semichord b.) The variables are used 
in this sense throu^out the rest of this paper. The kernel function 
(eq. (15)) can be written in terms of these nondlmenslonal variables as 


(15) 


ioacg 

K(xg,yo) = e ^ 
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K(xo^yo) = - 


2 


yi^c 


IkXr 


5 n(xo-piypi) ^. 

l/v - I 


_ 

|e u(xo-3|yo|)sin^^to2 _ p2y^2^ + 


k 

M, 


'p|yo! 


IkA 

e u(X-3|yo|)sln^^J?^T7%^jdA 


(l6a) 


An eqxil valent expression for K(xo,yo) which will he useful In subse- 
quent considerations Is 




IJ^kXg 

“P^/M-L.+ i 

\Mk Sxo M 


U (xq-P I yo t ) sin ^^Xq 2 - p2y^2 


k g-lkxo 
M 


pX . 

/p|°o( ® U(A-p|yo|)sln^^^ - p2y^2jdA 




(l6b) 


Isolation and Dlsciisslon of Singularities of Kernel Function 

As prevloiisly mentioned and as may be noted In equations (15) and 
(l6a), the kernel fimction becomes singular at yo = 0 and Is of an 
Indeterminate nature when Xq = ^|yQ|. It Is therefore desirable to 
make special treatment of the function in the neighborhood of these 
values of Xq and yo In order to be able to express the function In 
a form which Is more amenable for calcxilatlons . Ihe Indeterminate 
condition arises from the first term of equation (l6a) because of the 
manner in which the lanit f\anction has been defined for this analysis, 
(ihe denominator of the first term vanishes at Xq = Pyo- presence 

of the unit function in the numerator, however, renders this singularity 
indeterminate . ) 

In the next few sections the forms of the singularities are ex- 
tracted (see eq. (24)) and the aforementioned indeterminate forms of 
the kernel function are eJiplicltly determined (see eq. (20)). A form 
of the kernel function more suitable for calculation piirposes, since 
the troublesome points are isolated, is presented in equation (26) . A 
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mann er of integrating the kernel function across its singularities is 
given in equation ( 27 ) . Tbe singularities of the supersonic and subsonic 
cases are then compared. 

Indeterminate form .- Consideration is first given to the indeter- 
minate form, and it is convenient for this purpose to consider the value 
of K(xQ,yQ) at points on the positive branch of a hyperbola. (See 
sketch.) Ihe equation of the hyperbola may be written as 



In these equations e = 0 corresponds to Xq = pyQ, since elimination 
of 0 gives Xq 2 - p2yg2 = e2. 

After substitution of these expressions for Xq and Plyol into 
equation (l6a ) , the results may be written as 


K(c,e) - -sefs 


2--lke cosh 0 


u(ee-0)« 


Vl2g2 sinh^O 

_ ike cosh 0 

|u(€e-0)e 


ik0 cosh 0 


cosh 0 cos + 


cosh 0 


p2 siji ^ + 

IkA 


^ F e U(A-e slnh 0)sln(^l/A2 - sinh20^dA 

slnh 0 


(18) 

To obtain a limiting value of this equation for small values of e, the 
trigonometric and exponential terms can be replaced by terms up to the 
second power of e in a series expansion. If the result of performing 
this expansion is denoted by K’, the equation 
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K'(€,e) = 2 u(ea-6)( ' °°°‘* ° - 

V2^ leSsinh^e 


Ik 

P^€ 


k^cosh 9 _ k^ 2.og cosh 9+1 
2p^ 2p^ sinh 0 

Is obtained, which, in tenns of the original coordinates Xq and j-q, 
is 



K'(Xo,yo) 


Vl2 


u(xo-P|yo|) 


Xr 


\/xq 2 _ p2y^2 

k2iog isjLiSZ^S 

P|yc| 




2ik 



( 20 ) 


Althou^ these eqmtions were obtained in order to reveal the form of 
the indeterminate value of the kernel function, they are found to contain 
singxilarltles at yQ = 0. Prior to any further discussion of this result 
it is desirable to consider the limiting form of K(xQ,yQ) as yo 
approaches zero to determine all the singularities at yo = 0. 


Singularities at yo = 0.- For the pvurpose of obtaining a limiting 
value of the kernel function for vanishingly small valiies of - yo^ the 
Integral appearing in eqiiation (l6a) may be written as the sum of two 
integrals, namely: 



e u(A-p|yo|)sin^pJ?^T'PJ7jdA = 

_ ikA 

e U(A-p|yo|)sin^^l/?^~r’^^jdA - 

IkA 

U(A-P|yo|)sln^ppT72^dA 


( 21 ) 
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The first of these integrals may "be evaluated frctn the table of Laplace 
transforms of reference 17 and has the following value: 


ikA 

j ® U(A-P|yol)sln^^\f^"T^^jdA = M|yo[Ki(kjyo|) (22) 


where Ki is the modified Bessel fimction of the second kind, of first 
order. In the second Integral, the integrand may he replaced hy terms 
up to the second power of yg in a series expansion. Thus, 



ikA 

^ U(A-p|yo|) sin/pJ?Fri^jdX 



ikX 

e U(X-3|yo|) ^sin ^ cos 




( 25 ) 


where Ci and Si denote the cosine integral function and sine Integral 
function, respectively, which are defined as follows: 



ccfi.. dt 
t 


Sl(x) 




Bin t 
t 


dt 


Substituting equations (22) and ( 25 ) into equation (l6a) gives, as a 
limiting value of K(xQ,yo) , 
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11 m K(Xo,yo) 

yo~^o 


lim 
yo — >0 


f 


2 e-il^to 

VZ^o^ 


Hyo|Ki(i^|yol) 



Cl 



+ Si 




(xo-P|yo|) 



-(^) - -(S)]j 

where the following series expression for K^^Cz) (see ref. I 8 ) is 
employed; 


Ki(z) 


= (y + log + ^ + 

2 )\ 2 . 16 584-. 



+ 1 . + . 
z \4 64 1152 


(24) 



where y is Euler's constant (7 = 0.5772157)* Examination of equa- 
tion (24) shows that the only singular terms are the same as those which 

appear in equation ( 20 ) , namely - — _*.£ and -k^e” ° 

yo^ 

Thus, for the purpose of isolating the singularities of the kernel 
function, only K'(xQ,yo), as defined in equation ( 20 ), need he con- 
sidered. Nevertheless, the results given in eqiiation (24) may he useful 
in some applications since they provide a ready means for evaluating the 

nonsingular part of lim K(XQ,yo). 

yo — +0 



Form of kernel function suitable for calculations .- As in the suh- 
sonic case, with knowledge of the critical values of the kernel fvmction, 
an expression can he written in which the kernel function is separated 
into two parts, one of which contains no singularities or indeterminate 
values and the other of which contains all the singularities and critical 
values of the kernel function. This expression is 


K(Xo,yo) = K(xo,yo) - K»(xo,yo) + K'(Xo,yo) 


( 26 ) 
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where K(xQ,yo) is defined in equations ( 15 ) or ( 16 ), and K'(xQ,yo) is 
defined in equation (20). The term K(Xo#yo) - K*(xo,yo) in equa- 
tion ( 26 ) has no singular or indeterminate values. The term K’(xQ,yQ) 
is singular at yo = 0 snd Indeterminate when Xq = ^yQ. 

Integration of singularities of kernel fxmction .- Since integration 
of the kernel function is often necessary, a few remarks on how to cir- 
cumvent its Inherent singularities are in order. Each term of K’(Xo,yo) 
in equation (20) possesses a simple indefinite integral with respect to 
the variable T = y - yo- performing integrations with respect to t] 
that Involve a passage across the line H = y# a principal value is to 
be taken. For example. 



^(xo-3|yo[)xo dTi 

(y - n)^Vxo2 - p 2 (y _ ti)2 


xo 


U(xo-p|y-l|)/xo2 - 3g(y - 1)2 _ 


y - 1 


^(xo-P|y+l|) - p^(y + 1)^ 

y + 1 


( 27 ) 


where the symbol j' indicates that the singular integral is to be con- 
sidered simply as a function of its limits. A Justification for this 
consideration is that it leads to results that could, with considerable 
labor, be rigoroiisly established by maintaining the variable z in the 
analysis until all operations are performed. 


Comparison with singularities of subsonic case .- It may be of 
interest to compare the above results with corresponding results for 
the subsonic case, that is, for Xq > 0 and yp » 0. Results for the 
subsonic case may be obtained from eqtiation ( 5 I) ot reference 1 as 
follows : 


lim K'(xp,yp) = 11m ^ 

yo-^0 yp^O Yl^ 


^Xp2 + (1 - M^)yp2 + xp 

+ (1 - w^)yo^ 


ik 

l/xp2+ (l-I^)yp2 



k(]jxo^+ (1 - - Xp) 

2(1 - M) 


r__2_ + ^ Ipg ^^(3- + M)ypg 

Yl^ 2 4xp 


(28) 
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Hie singular terms of this expression for Xq > 0 and M < 1 are 




Comparison of this result with eqxiation (20) shows that the singularities 
for subsonic and supersonic flow are of identical form. 


Some Infinite-Series Expansions Pertinent to the Kernel Function 

The kernel function can he expressed as a series by various e^qiansion 
procedures. Some particular expansions, which should be iiseful in appli- 
cations, are discussed in succeeding paragraphs. These are the power- 
series expansion in terms of the reduced-frequency parameter (see eq. (29)) 
and an expansion in terms of Bessel fvinctions. The latter expansion is 
iised in a later section to obtain the Isemel fimction for two-dimensional 
flow from that for three-dimensional flow. 


Power-series expansion with respect to k . - As in the case of sub- 
sonic flow, the kernel function can be expanded into a power series with 
reswct to k that, in the present case, is useful for small values of 
k/p2^ a combination of reduced frequency and Mach nvunber that is Inherent 
in such an expansion of the supersonic kernel. The terms of the expansion 
may be simply obtained by e3q)andlng the terms of equation (l6a) that are 
fxmctlons of k and collecting the results. The first few terms are 


„ -ikxn 

K(xo,yo) =- ^%2 ^(^o~P|yol) 


^ Jtx 






- i)p%0^XQg -f (a - ^ 

slpsj |/xo2 - 


(29) 


Although this power-series expansion converges to the appropriate value 
of K(xo^yo) S'!! finite values of k/p2^ a great number of terms are 

required imless k/p2 Is small. These first few terms of the expansion 
can be considered to represent the kernel function for values of k in 
the range of magnitudes generally encountered in dynamic-stability studies 
and, therefore, they are pertinent for obtaining time -dependent stability 
derivatives. A noteworthy featiire of the expansion is that each term can 
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iDe Integrated, In the sense that it contains a simple indefinite integral., 
with respect to the variable tj = y - y©* When such integrations involve 
a passage across the line tj = y, a principal value is to he taken in the 
sense described after eq.mtion (28) . 


Expansions in terms of Bessel fimctions . - The trigonometric terms 
appearl^ in the expression for K(x^,yo) in equations (15) # (l 8 a), and 
(16b) can be expanded into infinite series Involving Bessel functions of 
the first kind. Such expansions have good convergence properties, even 
for large values of the parameter and each term possesses a simple 

Indefinite Integral with respect to tj. Such Bessel function series are 
therefore useful for deriving an expansion for the indefinite integral 
of K(xo,yo) with respect to i). The indefinite integral of K(xo,yo) 
leads to the downwash associated with pulsating vortex lines ("horseshoe" 
vortices) and, as previously indicated, to the kernel function for two- 
dimensional flow. It ml gh t he useful to point out that the expansion of 
the cosine term into a series involving Bessel functions is also useful 
for studying distributions of pulsating sources. 


For the purpose of expanding the trigonometric terms tmder dis- 
cussion, consider the expressions 


H(Va) 

- a 2 

and 

H(A-a) 

/a 2 . a2 

where a, b, and A are positive. 


U( A-a) s in ^b \/a^ - 


>)= 

^ 8 b 


(50) 


(51) 


By wialcing use of a known Fourier transform relation, expression (50) 
can be equated to an infinite Integral involving a Bessel function of the 
first kind (see, for exanple, p. 35 of ref. I7) : 



To(xf^) 


cos 


Jo 


(A^T^b^) cos 


ar dr = 


cos 


1Wa2 . S,s) 

- a2 


ar dr = 0 


\ 

(A > a) 


(52) 


(A < a) 
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By use of the addition formula for Bessel functions (see, for example, 
p. 558 of ref. 18 ), the Bessel functicm appearing in this equation can 
he written as an infinite sum of products of Bessel functions as follows; 

JoUf^) = Jo(tA) Jo(bA) + 2 ^ (-l)“j2n(T>^) J2n(^^) (55) 

' n=l 


Xhus, 


f Jo(x\(A2+h2) 

'^0 


cos aT dr = U(A-a) ■22S 


Uy? - »g) 


= r 




Jo('T>0 + 


n=l 


cos ar dr ( 5 ^) 


In view of the relation (see ref. IT, p. 5?) 

U(X-a)cos^2n sin"^ 


/ J2n('>•^) 

^0 


cos ar dr 




(55) 


the indicated integration on the ri^t-hand side of equation (jlf) can he 
carried out term hy term so that 


Jo(hX)+ 2 ^ (-I)“j2j^(h?0cosf2n sin"^ 
n=l ^ 

( 56 ) 

Substituting the expression on the right-hand side of equation ( 56 ) into 
equation (51) gives 


U(?v-a) cos(hVA^ - a^ j _ u( A-a) 

\/a2 - 


U(X-a)sln 




7«(7>-a) 

^2 - a2 

27Jj(X-a) 
_ a2 


^ oo8 ^2n sin"^ 

oos^Bin"^ |(2n - l)Bin ^ — J 


(57) 
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But, since 


U(A-a) ■ ■ ^ - cos fsln"^ = U(A-a) (a > 0) 


the expression for U(A-a)sin(hyA^ - a^j may he written as 

U(A-a)sin(b/A2 - a^) = 2U(A-a) ^ (-l)“"^J2n_i(l3A) cos 

n=l 

(58) 

By direct cco^jarison, equations ( 36 ) and ( 38 ) can he used to write 
expanded forms of the trigonometric terms appearing in equations ( 15 ), 
( 16 a), and (l 6 h) . Expansions thus obtained will now he used to derive 
the kernel function for two-dimensional supersonic flow. 

In contrast to three-dimensional flow, a physical interpretation 
of the kernel function for two-dimensional flow is that it represents 
the downwash at a given field point due to a pulsating hound vortex 
line of infinite length. IHiis kernel function may he obtained by inte- 
grating the kernel ftmction for three-dimensional flow from -w to 
or, in view of the role of the unit function, from one J&ch line to the 
other, with respect to the variable tl = y - yQ. Pulsating "horseshoe" 
vortices may he obtained by integrating K(xo^yo) over an arbitrarily 
finite range with respect to t^. 


( 


(2n - l)sln“^ - 
A 


Derivation of Kernel Function for Two-Dimensional Flow 

In this section the kernel function for three-dimensional flow is 
reduced to the function for two-dimensional flow, and the filial results 
of the reduction are given in equation (l}-9) • ^or the purpose of deriva- 
tion, K(xo,yo) will be considered as given in equation (l6b) . The 
two-dimensional kernel function can then be expressed as 



l^^l£Xo 


a +i 


poo 

j K(Xo,yo)<i’l = 

e ^ 

S Mk 

mI7 



T 

* «oo 

yo^ L 




= K(Xo) 


(39) 
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or 



In eqmtlons (4l) for and equation (42) for I 2 , infinite limits 

vlth respect to the integration of tj are not necessary. In view of 
the role of the unit function in the integrands, limits of +Xq include 
all values of q for which the Integrands are different frcm zero. Ohe 
symbol S indicates that the integrals are to he considered simply as 

functions of their limits or that the singularity at q = 0 is to he 
ignored. 

First consider equations (4l) for Iq and then perform an integra- 
tion hy parts. The expression for Iq may then he written as 
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By direct compaxison with eqmtlon (56) , equation (i|- 5 ) may be written in 
expanded form as 


2 ^ (-l)^J 2 n^p Xq^cob ^ 2 n sin"^ dr] ( 

In this equation, the terais involving. do not contribute to the 

Integral because 

r° eosfan XV, = 

■Jo . ,2 V xj 

i[^(xo-|«l|)sln(2.i sln-1 ° j^6(x„-|n|)sla(2n B±n^ i)dn 


= 0 


Hence, since 




n 


= I u(xo) 


the expression for can be written as 

. -X^ 8 (*o) + ^ °(*o) =ni.J - P ''(='0) Ji(p *0 




Now consider equation (k- 2 ) for I2; namely 

ikTv 




r"'°e''p2 u(A-|Ti|)slnf^^?^ Ti^jdA 

■^hl 


dTi 
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The double integral in this expression can be considered as a surface 
integral over a triangular region of the Aii-plane cut out by the 
lines T) = X, tj = -X, and X = Xq> as shown in the following sketch: 



By a change in the order of Integration, which is admissible since the 
singularity at ti = 0 is to be Ignored, the e:q>resslon for may be 

written as 

IkX , 

I 2 = P r ° e dX ^ slntej?^~I~7i2\dTi (1^6) 

^0 "'-x \p2 j 

The inner Integral in this eqmtlon is identical in fonn to the integral 
in equation (ij-l) . Hence, by observation of and comparison with the 
results obtained for in equations {k3 ) , (^) > and (45) f it is found 

that 

. -f: U(X)do(p x) (47) 

The ejipresslon for I 2 can therefore be written 

ikX 

Ig = ° e u(X)Jof— X^X (48) 

P 'Jq / 

Substituting this resiilt and the results given in equation (45) tor II 
into equation (40) gives a desired form of the Icemel function for two- 
dimensional supersonic flow: 
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K- (Xq) 



Examination of equation (it-9) shows that the only singularity involved in 
the kernel function for two-dimensional supersonic flow is the 6-function. 
At zero frequency, all the terms of K (Xq) except the 6-function vanish. 
Hie kernel function required to treat two-dimensional wings at steady 
angle-of -attack conditions is therefore proportional to this 6-function, 
and, as shown in the following section, leads in a very sinple manner to 
the well-known Ackeret results. 


Hie integral that remains to be evaluated in equation (49) is well 
behaved and similar to Integrals, treated by Schwarz (ref. I 9 ) and others, 
that arise in the velocity-potential approach for treating two-dimensional 
wings . 


Application of Kernel Function to Lift Distributions 
for !I!wo-Dlmensional Wings 

Hie results obtained in the previous section for the two-dimensional 
kernel function are now onployed to obtain the lift distribution on 
oscillating and steady two-dimensional wings moving at supersonic speed. 
(See eqs. (56) and (61).) Since the lift distributions so obtained agree 

with the Ackeret results for a steady wing and also with known results 

for the oscillating wing (ref. 15 ) ^ they serve as a check on the correct- 
ness of the expressions for both the two-dimensional and three-dimensional 
kernel functions. 

Hie integral equation that must be solved to obtain the lift distri- 
bution for two-dimensional wings in supersonic flow is particularly 
simple since it involves a single integral of the convolution type; 

' tfe r jC ‘50) 

Integral equations of this type can be readily solved by laplace 
transform procedures since the laplace transform of a convolution integral 
is the product of the transforms of the functions that ccmpose the 
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IntegrsLnd. In the present case. If s represents the laplace transform 
operator defined by 


L|f(x)J = J f(x)dx = f(s) 

the transform of equation ( 50 ) may be written as 


(51) 


w(s) = ^ L(s)fc (s) 
4itp 


(52) 


Solving this equation for L(s) gives the Laplace transform of the lift 
distribution: 


L(s) = 

I 1C (s) 


(53) 


Inversion of the transform on the ri^t-hand side of this equation gives 
the lift distribution- 

For the case of a steady two-dimensional wing. 


w(x) = Va 

w(s) = ^ 

K(x) = ^ 5(x) 
VI 

K(b) 

VI 


Then 


Us) 

^s 


(54) 


The inverse transform of equation (54) gives for the lift distribution: 

(55) 


L(x) = U(x) 

P 


From this result, the total lift per unit of span is 

Chord 

L(x)dx = 

P 


P Chord 

2/ 2 L(x)dx = 

^0 B 


(56) 


This result agrees with the well-known Ackeret result. 
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New consider the unsteady case for oscillatory translation, 

w(x) = lojlh = iVlsh (57) 

where E is the ampUtiide of displacement referred to l, and the 
Laplace transform of w(x) Is 


w(s) 


IVkfa 

s 


(58) 


The laplace transforms of the different terms of /C(x) (eq. (^ 9 )) can 
he simply derived or they may he obtained from Laplace transform tables 
(for example, ref. I6) . After conODlnlng the transforms of the different 
terms, the restilts can he written as 


/C(s) 



M2k2 


V 2 


s + Ik 


( 59 ) 


Substituting equations (58) and ( 59 ) into equation ( 55 ) gives for the 
transform of the lift distribution 


L(s) = 8 + ik 


s /s + 


ii^k 


^ l^k 2 


(60) 


The Inverse of this transform gives for the lift distribution 


L(x) = aievM 


U(x)< 


it^kx 

P 2 


’iW) • “X" 


il#kg 

U(|)e' 




(61) 

This resiilt can easily be shown to check with the resiilts of refer- 
ence 15. Moreover, if il£ is set equal to a, and then k is allowed 
to approach zero, equation (61) reduces to the res;ilt for the steady 
case. 


CONCLUDING REMAEKS 


The main purpose of this paper was to derive and present in a form 
that could be nimierlcally evaluated the kernel function of the integral 
equation relating downwash and lift distributions for oscillating wings 
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in sijpersonic flow. This purpose has heen achieved for •bhree-dimensional 
flow, and the results have been converted to a fona more sviitable for 
calculation by isolating the singiilar or critical points. The kernel 
function for two-dlmenslonai supersonic flow has been presented and the 
results show that the only singularity is a Dirac delta function, which 
appeared in such a manner that further reduction with regard to slngu- 
lEiritles is not required. 

The results presented in this paper for supersonic flew together 
with those previously obtained for subsonic flow provide a kernel function 
that is capable of being evalmted at any Mach number. As experience 
develops it is expected that use can be made of the kernel function to 
develop approximate procedures, that will be more or less uniform throu^- 
out the Mach number range, for calculating aerodynamic forces on oscil- 
lating (or steady) wings of arbitrary plan form and with arbitrary down- 
wash conditions. The labor Involved in such approximate or numerical 
procedures will Indeed be prodigiotis and will require the use of modern 
high-speed computing equipment. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
langley Field, Va. , Pebriiary I 5 , 1955* 
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APPENDIX A 

REDUCTION OF KERNEL FUNCTION TO SONIC CASE 


The purpose of this appendix is to reduce the results obtained for 
the kernel function in sijpersonic flow to the sonic-flow case. Such a 
reduction shoxild, hy ccmparison with results obtained for the sonic case 
in reference 1, provide a partial check on the correctness of the re^^ts 
for the supersonic case. To effect this reduction, it is first necessary 
to discard terms arising from the advanced-type potentials employed in 
deriving the doublet potentials for supersonic flow. This may be accom- 


, J^Jxo2-32yo2 

plished by replacing cosf^^VxQ^ - with ^ 


and 


\3' 

i^l/x^ 2 _ 32 y 2 


0 


in equation (l6a) of the 


analysis. The limiting value of K(xQ,yo) as M approaches 1 can 
then be written as 


K(Xo,yo)j^l 


= lim - 


^-ikxo 

XoU(xo-Plyol) 

Vl^o^ 

i/xq 2 - 32y^2 


Ik 


^Xo-M|/xQ2_32y^2 ^ 


I U(x^-3lyol ) 


^(xo-M||xo2_p2y^2j 


pxo 


(Al) 


When the limit M = 1 is approached frcm the supersonic side, the term 
M is conveniently replaced by 


M = 1 + — e 


(A2) 


where e is infinitesimally small, so that 

32 = (M - 1)(M+ 1) = |e ^2 + | 6 ^ 


as € 
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With this approximation, equation (A1) can he written as 


-ikxo 

K(Xo,yo) = lim 

e— ^0 VZ^yo^ 


XoU(xQ-,/elyol) 

|/xq 2 - ey^2 


-^U(xo-€|yol)e“ 


/Xo^-eyo^ 


ik 


1 + l^/iyo 




€ 




= Iim _a 


-ikXo 


e— ^0 VZ^Yo^ 


2U(Xo)e 




V 

nxo -—[a 

-fl+£) 

k- . 

/ U(A)e ^ L 

\ 2/ 

\ 2A Ij 



* 


VZ^y 2 


o i_ 


2U(xo)e 




ik 


r^o 

/ U(A). 

OQ 


■ ^o' 

- -y. 


dA 


(A5) 

A next step is to show that the result in equation (A5) is equivalent 
to that given in reference 1. For this purpose, the final integral in 
equation (A5) is considered as two integrals: 



(A4) 


The second Integral on the ri^t-hand side of eqiiation (a 4) is already 
in suitable form. The first integral on the right must he converted as 
follows : 
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Since the argument of the unit function in the integral on 

the rigjht in equation (A5) is always greater than zero for yo ^ 0, this 
integral has the same value as the integral 



Integrating this expression hy parts with 



(A6) 

From, reference I 8 (pars. 6 .I 6 and 10.41), this result may be written as 



52 


NACA IN 3458 





+ Yq 


= i 

2/ ^ 


-i + iyo' 


L 


cos 


0 (52 + y^2)5/2 


+ 



sin 


^ 1 
k 


-i + ik|yo|Ki(k|yo|) - 


|i'|yo|[ii(!=|yo|) - ii(''|yo|) 



where Is a modified Bessel function of the second kind, I-j_ is a 

Bessel f\jnction of the first kind with an imaginary argument, and 
Li = lui - ^ is a Struve function of the first kind with an imaginary 

argument. Substitution of the second integral in equation (a 4) and the 
results of equation (A?) into equation (A3) gives, for the sonic case. 


K(Xo,yo)j^l " 2 j -2U(0Co)e 


2 




+ 1 + ikyo - 


klyo|Ki(k|yol) - ^ |yo| [li(k|yo|) - Ii(k|yo|)J 


ik 


f ° U(A) 

^|yoi 


ik A yo2\ 


(A8) 


A cojnparison of this result with the result given in equation (47a) of 
reference 1 shows that the two equations are equivalent. 
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APPENDIX B 

DERIVATION OF DOWNWASH FUNCTIONS ASSOCIATED WITH 
"HORSESHOE" VORTICES IN SUPERSONIC FLOW 


T^<a dcwnwasb. associated with a vortex line can be obtained by an 
Integration, between appropriate limits, of the kernel function K(xo,yo) 
with respect to tj = y - yg. In order to perform such an Integration 
analytically, recourse must be had to term-by-term Integrations of an 
expanded form of K(xQ,yo) . In this regard, use can be made of the 
expansions given In equations (38) and (38) of the analysis to obtain 
expanded forms of the downwash functions for vortex lines that have veiy 
good convergence properties, especially for the range of values of the 
parameter Mk/p 2 that would usually be of Interest In applications. 
Expressions so obtained will be cumbersome and will require high-speed 
conputlng equlEment to make them very \iseful. 

Tn regard to "horseshoe" vortices In supersonic flow, there are 
five different significant regions In which a field point may be con- 
sidered to be located (see sketch) . 


1 ^ 




y 

y 

w 

y 

Trailing vortex 


( 5 ) 


\ 


X 

§ 

t» 


\ 


\ 


( 2 ) 












\ ' Wind direction 

(1) > (3) »- 

/ \ 




\ 


Trailing vortex 
\ 


( 5 ) 
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Region (l) is between tbe Mach cones emanating from the end points 
of the bound-vortex line. The downwash at a point in this region is not 
affected by the trailing vortices but is created by the bovind vortex 
alone. Therefore, the downwash is the same as would be produced by a 
bound vortex of infinite length and corresponds to the kernel function 
for two-dimensional flow discussed in the analysis. Region (2) is 
between the trailing-vortex lines and is within the Mach cone emanating 
from one end of the bound vortex but outside the Mach cone emanating 
from the other end. !lhe downwash at a point in this region is created 
by the bound vortex and one of the trailing-vortex lines. The other 
trailing vortex has no effect on the downwash. Region ( 5 ) is between 
the trailing-vortex lines and is within the Mach cones emanating from 
both ends of the bound vortex. Downwash in this region is created by 
the bound vortex and both trailing-vortex lines. Region (4) is outside 
the trailing-vortex lines and is within the Mach cone emanating from one 
end of the bound vortex. The downwash is created by the bound vortex 
and only one of the trailing-vortex lines. Region ( 5 ) is outside the 
trailing-vortex line but within the Mach cones emanating frcm both ends 
of the bound vortex. The downwash is created by the bound vortex and 
both trailing-vortex lines. 

For any of the five regions discussed in the preceding paragraph, 
the integral corresponding to the downwash function may be formally 
written, with use of equation ( l6b) , as 



where use of the substitution Pyg = ^ gives 
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and 


3(y-ii2) 


I 


Kl 


u(A-|^I)sln^pl(F 


^2]dA 


(B5) 


In eq\iations (B2) and (B5)., a principal part - as described after eqim- 
tion ( 27 ) in the analysis - is to he taken when the integrations are 
carried across the line ^ = 0- Ihe purpose now is to reduce these 
expressions to forms amenable to nvunerical evaluation. The first step 
in this procedijre is a reduction of the expression for % (eq. (B5)) • 
The double Integral in this expression can be considered as a surface 
integral in the -plane where the order of integration is first with 
respect to A and then with respect to t,. The steps in the reduction 
are first to delineate the area of integration for each of the five 
different cases imder consideration, and then to change the order of 
Integration in the surface-integral representation of Ijj.. 

From the description of the different cases to be considered and 
by examination of the limits of integration in equations (B1) , (B2) , 

(B3) , the area of integration for ilie case of a field point in each 
of the aforementioned regions may be considered as shown by the hatched 
areas in the following sketches; 
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Expressions for Ijj. for the five different regions or cases may then 
he e^qjressed as sin 5 >le integrals as follows: 


Case ( 1) : 


Case (2) : 



ikA 

e F(A,-A)dX 



IkA 

F(A,-A)dA 


X 


^1 


ikA 

P^ F(A,-5l)dA 


(B4-) 


(B5) 


Case (5) : 




ikA 

e"p2 


F(A,- 



F(X,-^l)dA 



_ik2i 


Case (4) : 



IkA 

F(A,^l)dA 


(B6) 

(B7) 


Case (5) : 



ikA 


F(A,^l)dA + 


X 


^2 


_ 


(B8) 


The expression 

F(A,-A) = ) sin ^^a2 - ^2 


is evaluated in the text in ccmnecticai with the derivative of the two- 
dimensional kernel function and is found to reduce (see eq. (47)) to 
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F(-K,-X) = ^ U( X) Jq a) (bio) 

The F-fxmctlon for other arguments can he obtained hy substituting appro- 
priate limits in an integration by parts of namely (see the 

development following eq.. (li-l) in the analysis) : 

F(X,-^l) sin p - ^2 d^ 



After the first term on the ri^t-hand side of equation (BU) has 
been expanded by coo^arison with the expansion given in equation (58) of 
the analysis and the limits of integration have been substituted, this 
expression may be written as 






m 



00 



n=l 


2(-l) 


n-1 


5l 


T- /MkA\ Vf 


(2n - l)sin‘ 


-1 


3^ n 




+ 


(B 12 ) 
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Substitution of ^2 ^ ^ limits of equation (Bll) gives 

F(^2,-^i) = -u(A-|?2|)j 


m j te\ sin-1 + 

p2 0i^p2 j X 


n=l 


2(-l) 


n-1 


[(2n - Dsin"! ^ 


Mk (-1)^ 




2 n r2 


(m) 


sin|2n sin 


-Ik' 


} \ 


u(x-l?i|) 
2(-l) 




n=l 


p2 ; 

n-1 


^2n-l 


/MkX^ 

vp^y 


cos 


(2n - l)sin“ 


'9 


2n sin-1 k 


Mk (-1)“ T (mA (, 

If -^1 1® replaced by limits of equation (BU) , then 


F(X,?1) = 


jrMk 

P2 


n=l 


U(X)Jo(S2^) . U(v|y){p ^ 


j, 




(BI 5 ) 


(Bl 4 ) 


Substitution of ^2 


for 


A and for in equation (Bll) gives 
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When equation (BIO) equations (B 12 ) to (BI5) are suhstltuted into 
equations (B 4 ) to (B 8 ) , respectively, they give the reduced forms of fi,' 
lor the five cases xinder consideration. ” 

After the reduction of the corresponding reduction of Ix is 

considered. As may he found hy examination of the eaqu-ession for 
(eq. (B 2 )) and the sketches showing the areas of integration for the 
^i^^’srent cases, reductions of corresponding to those of Ijj, can 

he obtained fran the P-functions (eqs, (BIO) to (BI5)). Results for the 
different cases may he expressed as follows: 

Case (1) : 

Case (2) : 



4o 
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Case ( 3 ) : 


Case (4) : 


Case ( 5 ) : 



(BI 8 ) 


(BI 9 ) 


(B20) 


When the expressions for (eots. (B16) to (B20)) and 

(eq,s. (b 4) to (b 8 )) that are associated with each particular case are 
substituted into equation (B1) , expressions for the downwash at each of 
the five significant field-point locations may he obtained in terms of 
the F-functlons (eqs. (BIO) to (BI 5 )) as follows: 

Case ( 1 ) : 


K(Xo,yo)dT = e 11m ^ 


il^kxr 


A — >Xq 


k -itoco r°.“^ 


M 


Case ( 2 ) : 


r K(xo,yo)dT] = 

^Til I 


^0 


f il^kxo 

p— 


ihA 


F(A,-A)dA 


(B21) 


lim 


k g-ilKo 
M 


A— >Xo\Mk 8A 

~ t - 

f ^ e F(A,-A)dA + 

•Jci 




r^o “ q2 

/ e P F(A,-^l)dA 


0 
IkA 


(B22) 
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Case ( 5 ) : 


f ^ K(Xo,yo)dn = 


k e^ikxo 

M 


r l#kXQ 

Um ^ ^ 

A — >Xq 

/ e 3^ F(X,-A)dA-f 
Jci 


-V(?2 -^1) + 

\Mk: 8X My ^ 


'^1 


IkX 

e 3^ F(A,-^ 3 ^)dX + 


f, XX.A 


ikA 


'^2 

Case ( 14 -) : 

f ^ K(Xo^yo)<3-n = 


r o o2 

J e ^ F(C2,-?l)<i>^ 


(B25) 


ll^kx p 

I 32 lim 


X »-Xr 


IkX 




M 


J" ° e 32 F(A,^i)dx[ 


Case ( 5 ) : 
'^2 


f ^ K(Xo,yo)<iT] = -^ 

^7,1 Z 


5l 

lly^kxo 

Um /JBS ^ + 1 ' 
X >Xq 


(B214-) 




ke-liX; 

M 


L 


^2 - — 

^ e F(X,5l)dX + 

^1 


C 


IkA 


F(C2,^l)<iZv 


(B25) 


The results fdr case (l) (eq,. (B21)) agree with results obtained for the 
two-dimensional kernel function given in equation (14-9) of the analysis. 
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